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Abstract 

, In this paper, we are first interested in the compressible Navier-Stokes equations with density- 

' dependent viscosities in bounded domains with non-homogeneous Dirichlet conditions. We study 

, the wellposedness of such models with non-constant coefficients in non-stationary and stationary 

(-H ' cases. We apply the last result in thin domains context, justifying the compressible Reynolds 

^ , equations. 
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Introduction 



The Reynolds equation is a linear equation describing the evolution of the pressure in mechanisms of 
lubrication. More precisly, it is used to calculate the pressure distribution in a thin layer of lubricant 
between two surfaces. It was proposed by O. Reynolds in 1886, see [23]. This equation is very much 
used in mechanics, for instance to describe the process of lubrication of magnetic hard discs. In 
] the same way, air flow between the two surfaces constituting a rigid disk assembly (flying head and 
' magnetic storage surface) is frequently modeled using the Reynolds equation. 

It was proven, first in 1986 by G. Bayada and M. Chambat, see [2], that the Reynolds equation is an 
approximation of the Stokes equations in thin cases. This proof was formulated by taking as initial 
Q\ ] equations the incompressible Stokes model. Since, many other works (see |21] and the cited references) 
has made it possible to refine the first result by giving errors of the approximations between the Stokes 
t> ' model and the Reynolds incompressible model. 

^ ' In the previous examples of applications the fluids (like air) are clearly compressible fluids. Within 
5^ the framework of the compressible fluids, there exist a so called compressible Reynolds equation which 
is, at least formally, the asymptotic of the Navier-Stokes compressible equations in a thin domain. 
Contrary to the incompressible classical Reynolds equation, the compressible Reynolds equation is 
highly nonlinear and has been a subject of many mechanical studies [3l |12] or of numerical studies [H 

m- 

However the literature about the rigorous justification of these equations in the compressible case is not 
very important. It would seem that there is only one result, due to E. Marusic-Paloka and M. Starcevic 
(see [18] or more recently [19]), restricted to the case of ideal gases. The primary reason of this lack 
of literature certainly comes from the fact that the study of the compressible Navier-Stokes equations 
is rather difficult. Recent works of D. Bresch and B. Desjardins on these compressible Navier-Stokes 
equations, see [7] for instance, showed that there exists a particular structure to these equations. 
In this article, we adapt these new results to use them and rigorously justify the compressible Reynolds 
equation for rather general state laws. 
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More precisely, we prove the two following results (precise statements are respectively given on page [5] 
and page [T7|) : 

Theorem 0.1 There exists a steady-state solution to the compressible Navier-Stokes equations in a 
bounded domain, with Dirichlet boundary conditions. 

Theorem 0.2 The compressible Reynolds equation is an approximation of the stationary compressible 
Navier-Stokes equations. 

The present paper is composed of the following parts: 

• In the first section, we present the notations and the classical compressible Navier-Stokes equa- 
tions. We give the assumptions as well as the theorems related to the compressible Navier-Stokes 
equations (in the non-stationary case and in the stationary case). 

• Section [2] is devoted to the proof of the existence result for the compressible Navier-Stokes 
equations in the non stationary case. 

• Section [3] is devoted to the proof of the existence result for the compressible Navier-Stokes 
equations in the stationary case. 

• In Section m we introduce the lubrication problem in term of thin film flow. We also anounce the 
convergence result for the stationary compressible Navier-Stokes equations to the compressible 
Reynolds equation. 

• Section [5] is devoted to the proof of this convergence result. 

1 Wellposedness of compressible Navier-Stokes equations with 
density-dependent viscosities 

1.1 Statement of the problem 

Compressible Navier-Stokes equations: The compressible Navier-Stokes equations describe the 
evolution of a compressible fluid in a physical domain i7 C M.'^, d G {2,3}, via the conservation 
equations of the mass and the momentum. They thus couple the velocity u of the fluid and its 
density p: 

dtp + div(pu) = 0, 
dt{pu) + div(/7U (g) u) = div((T) + f. 

To close this system, we must give the force term f and the stress tensor cr. 

• The force term f allows to represent friction forces or corresponds to a turbulent drag force. 
They read 

f = -rop\u\u, 
where rg is a non negative real coefficient. 

• Finally, we give the rheological law for the stress tensor cr : the fluid is assumed to be Newtonian, 
so that there exists two viscosity coefficients (called Lame coefficients) p = p{p) and X = X{p) 
such that 

a = 2pD{u) + (Adiv(u) - p)ld. 

In this equation, D{u) corresponds to the strain tensor (the symmetric part of the velocity 
gradient tensor), and the pressure p is determined using a thermodynamic closure law, i.e. an 
explicit relation p = p{p). 
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The equations in which we will be interested are thus: 

dtp + div(pu) = 0, (1) 

dt{pvL) + div(/)u u) + Vp{p) = div(2^(p)D(u)) + V(A(/5)div(u)) - ro/5|u|u, (2) 
and also the stationary correponding ones: 

div(pu) = 0, (3) 

div(/9u u) + Vp(p) = d\\{2p{p)D{\y)) + V(A(/9)div(u)) - rop|u|u. (4) 



Boundary conditions: The physical boundary conditions which interest us here (see part [4]) are of 
the nonhomogeneous Dirichlet type on the velocity field. For technical reasons, we will also need to 
impose a condition on the density on the boundary. The conditions will thus be the following ones: 

u = Ufo a given function on dVt such that • n = 0, 
p = Ph constant on each connected component of 50. 

The biggest part of the published works concerns the whole space case = R^, or the periodic case 
$7 = (see for instance [3 [9]). More recently in [8], the authors deals with the Dirichlet homogeneous 
condition or Navier's condition on the velocity field. The building that we present here draws hard 
inspiration from this last paper. Particularly, the boundary condition on p, being already present 
in [8], it is not amazing to find it in a more general case. 

Initial conditions: In the non-stationary case, it is necessary to give the initial conditions corre- 
sponding to the situation at time i = 0. The physical quantities for which we give information are the 
density and the momentum: 

p\t=o = Po and pu\t=o = uiq. (6) 

We note that, due to the non-penetration condition u;, • n = 0, integrating with respect to the spatial 
variable the mass conservation equation ([T|) we obtain 

Consequently, the quantity p does not depend on time and will be denoted Mq . 

This system ©-(H]) has been widely studied, starting from the case of constant coefficients A, p and 
pressure laws of type p{p) = ap^ (see notably [13l HH [151 US [12]). More recently, many studies have 
focused on density dependent viscosity coefficients A = A(p), p = p{p) in space dimensions 2 or 3. 
These studies were originally developed on Korteweg and shallow water models, corresponding to 
7 = 2, A(p) = and p{p) = p, see [H [51 [TJ [51 12] • They all rely on a new mathematical entropy (the BD 
entropy), that has been discovered in its general form in [7]. It requires that the following algebraic 
relation holds: 

Vs > 0, \{s) = 2{sp{s) - p{s)). 

We introduce in the next part, the hypotheses which we shall use later. Obviously, these hypotheses 
take back principally those of papers named here. 



1.2 Assumptions 

Concerning the viscosity coefficients A and p, we assume that A and p are respectively C'^(M+) 
and C^(R+) and satisfy 

Vs > 0, \{s) = 2{sp'{s) - p{s)). (7) 
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We also suppose that /i(0) = 0, that there exists positive constants cq, ci, Cq, (^i,A,m> 1 and | < n < 1 
such that 

Vs G]0, A[, CO < fi{s) < -s", c'q < i^'is) < (8) 

Co Cq 

Vs £]A, +oo[, cis"^ < < -s", c'lS™-^ < /x'(s) < \s"'~\ (9) 
Finally, we are interested in a pressure term of the following form 

P{p) = Ph{p) +Pc{p), (10) 

where Ph{p) = o-p'^ (a > and 7 > 1) corresponds to the classical equation of state whereas Pc{p) is a 
"cold" component. We assume that there exists positive constants C2, C3, p^^, (5 and q > 1 such that 

Vpe]0,p4, ^^^^'c(p)<;^, (11) 

VpGK,+oo[, -^2:|^<p^(p)<c3P^-^ (12) 

Recall that such assumptions were initially introduced in [7J in the framework of barotropic flows. 
Moreover, in the three-dimensional case, we impose (in fact in the stationary case result) 

m<7 + n-^, /3<2(7 + n-l), (13) 
and (to control the "thin domain" -dependency) 

7 

m < a — n + -. (14) 



1.3 Existence results 

Definition 1.1 We shall say that (/O, u) is a weak solution of ([i])-(0j with boundary conditions (0) if 
it satisfies following regularity properties 

peL^{0,T;L^n)), VpueL~(0,r;L2(Q)), ^ip{p) e L'^{0,T; L''{n)), 



VMpJvu g L2((o,r) X n), pu^ G ^^(0,2") X n), 

where ^ being taken such that ^(/o) = p for p < p^/2 and £^{p) = for p > p^,, as well as boundary 
Dirichlet conditions on u in L'^{0,T; L^{dQ)), boundary conditions on p in L^{0,T; L°°{dQ)), and 
equations (QP-fD in V'{{0,T) x n) for all T > 0. 

As usual, we deduce from these regularities and the Navier-Stokes system itself that p and u are 
continuous in time with values in W~^'^{Q), which allows to define their initial values. 

Theorem 1.2 (Non-stationary case) Assume that conditions ^)-lil2\) are satisfied and consider 
some functions po ^^^^^^ n^o such that 

Po Po 

where xQ"{x) ■.= p{x). 

Then, for all ro G there exists a weak solution of the system 

dtp + div(pu) = 0, 
dtipu) + div(pu u) + Vp{p) = div{2p{p)D{u)) + V(A(p)div(u)) - rop\u\u, 

associated with the initial conditions in the sense of the Definition \l. 1[ 
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Definition 1.3 We shall say that (p, u) is a weak solution of with boundary conditions if 

it satisfies following regularity properties 

where ^ being taken such that ^{p) = p for p < /5*/2 and ^(p) =0 for p > p*, as well as boundary 
Dirichlet conditions on u in L^{d^l), boundary conditions on p in L°°{dU), and equations 
in V'{n). 

Theorem 1.4 (Stationary case) Assume that conditions ^}- il3\} are satisfied. 
Then, for all rg G M^, there exists a weak solution of the stationary equations 

div(/9u) = 0, 

div{pu(g)u) + Vp{p) = div(2p(/))L>(u)) + V(A(p)div(u)) -rop|u|u, 
in the sense of the Definition \1.3[ 

2 Proof of Theorem 11.21 (non-stationary case) 

The main idea is to obtain good energy estimates, notably using the BD entropy, structure discov- 
ered by D. Bresh and B. Desjardins in [4]. This will provide enough compactness on a sequence of 
approximate solutions to pass to the limit and obtain a global weak solution. 

More precisely, the first step is to obtain suitable a priori bounds on {p,u), and next to consider 
sequences (pk,Uk) of uniformly bounded weak solutions constructed from an adapted approximation 
process. Such sequences may be built by using the regularization scheme given in Section 12.31 (see 
also ^). It leads to regular approximate solutions, still preserving physical bounds and the mathe- 
matical entropy, uniformly with respect to smoothing parameters. 

The scheme of the proof will be therefore the following. In Section 12.11 we recall the main idea of 
the Bresch-Desjardins strategy, we then deduce energy estimate and so called BD estimate (Subsec- 
tion [22]) • In Subsection 12. 3t we give the construction of approximate solutions. 

2.1 Bresch-Desjardins strategy 

The BD entropy is the dedicated idea to get many wellposedness of non-stationary models with non- 
constant coefficients, for instance the compressible Navier-Stokes equations (see 0) and it has been 
recently enlarged to some neighbour contexts like Shallow- Water (see [6]) or Magnetohydrodynamics 
(see i24j). 

The particular point of this strategy is the mixture between the mass equation and the momentum 
equation to get the control of non-linear diffusive terms with density-dependent coefficients through 
the BD entropy. 

We first multiply ([I]) by ip'{p) = ^^-^ to get: 

P 

dtifip) + u • Vifip) + /i'(p)div(u) = 0. 
Then we derive with respect to the space variables: 

dtVfip) + (u • V)V^{p) + Vu : Vfip) + V(/i'(p)div(u)) = 0. 
Let's now multiply by 2p, then, noting U = 2V(p{p) and using ([1]), we write: 

dtipV) + div(/>u U) + 2Vu : Vp{p) + 2/)V(/i'(p)div(u)) = 0. 
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Rewriting the last two terms, one has 

2Vu : V/i(/j) = 2div(/u(p)Vu) - 2;u(p)V(div(u)) 

= 2div{n{p)D{u)) + 2div{p{p)A{u)) - 2V(/i(/9)div(u)) + 2V/i(/9)div(u), 
2pV(/i'(/9)div(u)) = 2V(/9^'(/9)div(u)) - 2V^(p)div(u), 

where we recah that D is the symmetric part of the gradient, and where A is the skew symmetric part 
of the gradient. Then, summing with the momentum equation ([2]), we get 

dt{p{u + U)) + div(/)u (u + U)) + Vp{p) = div(2/i(/9)^(u)) - rop\u\u 

+ V{{2pp'{p) - 2/i(p) - A(p))div(u)). 

Notice that the assumption ([7]) is now necessary to get the following interesting form: 

dt{p{u + U)) + div(pu (u + U)) + Vp{p) = div(2^(/>)^(u)) - rop|u|u. (15) 

2.2 A priori estimates 

To control the boundary terms in various integrations by parts, we introduce a lift of the velocity field. 
Let u be a regular function such that 



u = u on u • n = on 50 and div(u) = on Q. 



2.2.1 Energy 



The energy estimate comes from the multiplication of ([2]) by u — u. Using equation ([T]) and the 
boundary conditions on u and u we obtain 

^ [ (p\^ + Q{p))+ [ 2f^{p)\D{u)\' + [ A(p)|div(u)p+ro / p\uf 

"■1^ Jn \ ^ J Jn Jn Jn ^^^g^ 

= ^ / (/)u-5)- / (pu®u):V5+ / 2p{p)D{vL) : D{u) + TQ [ p\u\u-u, 

It Jq Jq Jq Jn 

where Q'{p) := n(p) and pn'(p) := p{p). 
2.2.2 BD entropy 

The BD entropy estimate comes from the multiplication of (|15p by u — u + U. Using equation ([T]) and 
the boundary conditions on u, u and p we obtain 

jt i {p^-^^-Y^ + ^(^)) + X 2/^(^)1^^")'' + X ^^^^^ ■ ^^^^^ ^ '° I ^ 



dt 



[ {piu + U) • 5) - [ {pu0{u + U)) : VS - / 2fiip)A{u) : A{u) (17) 
Jn Jn Jn 

+ ro p\u\u-u-ro / /9|u|u • U. 
Jn Jn 

One of the main interests of this estimate is that not only it makes it possible to have a control on U 
via the control of p{u + U) but also that the "pressure" term Vp{p) ■ Vip{p) is very rich. 
Separating the pressure into two terms : p = Ph + Pc see assumption pU]) . we write 

/ Vp • Vifip) = I vphip) ■ Vifip) + / vpcip) ■ y^{p) =■■ h + Ic 

Jn Jn Jn 
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About the term I^-, we use the definition oi ph{p) = ap^ and find 

h = aj [ ^,'{p)p^-'\Vp\^ 
Jn 

We write the term Ic as follows 

Jn P Ja P 

Using the assumptions (jlip and (jl2p we obtain 



/c>- / p-"-V(p)|vH^W-^ / p^-V(p)|Vp|^. 

Moreover, for p < A we can use the assumption dS]) on p' and deduce 



Adding the two contributions Ih and /c we obtain 

where M := -'^-'^+^ and ^ being taken such that ^(p) = p for p < ^min(p*,74) and ^(p) = for 
p > min(p^,, A). 

Note that the contribution ^ p'^~'^ n' {p)\'V p\'^ allows us to control positive power of p. In fact, using 
assumptions ([5]) and that is separating small and large densities, we show that 

/ p^-V'(p)|Vpp >Co [ p-'-'+^lvpi^ 
Jn Jn 

where Co depends on Cg, c'^ and We finaly obtain the following inequality. 



(18) 



where M := -°-/+" < and iV := > 0. 



2.2.3 Control of integral terms 



In this part, we are particularily interested in the case tq = 0. In this case, the evolution terms are 
enough to control all the other terms. Nevertheless, an additional friction term with ro > naturally 
preserves the following calculations. 

Integrating with respect to the time t G [0, T] the energy estimate ()16p . we obtain 

|u(r)P ^ rT r 



p{T)^^^ + Q{p{T))) + [ [ 2p(p)|Z)(u)p+ [ I A(p)|div(u)|2 

/ Jo Jn Jo Jn 

+ Qipo)) + / p(T)u(T) • S - / mo • S (19) 
/ Jn Jn 



I mop 
2po 



[ /(pu®u):V5+/" [ 2fi{p)D{u) : D{u). 
Jo Jn Jo Jn 



7 



Each term of the right hand side member is controlled as follow: 

|u(T)|2 



p(T)u(T) • 5 



^ ^ Jn 



— n \- moo^m 



I / (pu (g) u) : Vu 



< IV5|, 



T 







T 



Jn 



2/i(p)Z)(u) : D(G) 



< 



< 



T 



Jn 

T 



Mp)|I?(u)P+ / / M/^)|I?(5)| 



Jn 



JO 



Mp)|I?(u)|2 + |Z)(5)|L / U{p) 



Jn 



We deduce from (fT9|) that 
|u(T)p 



+ Q(/,(T))) + / / /x(p)|I)(u) 
/ JO Jn 



Jn 



A(/^)|div(u)|^ 



<C(mo,po,u) + |VuU / / Hu|' + |I?(u)|L / / ii{p). 

Jo Jn Jo Jn 



(20) 



In the same way, integrating with respect to the time t G [0,T] the BD entropy estimate p7|) . we 
obtain (recall that in this subsection the friction coefficient is assume to be zero) 



Jn 



+ / / 2Mp)|^(u)|2+ / / Vp{p)-V^ip) 



Jn 



Po- 



|uo + Uol 



+ QiPo)) + [ (p(T)(u(r) + U(r)) • 5) - / (po(uo + Uo) • 5) (21) 
/ Jn Jn 

- / / (pu ® (u + U)) : V5 - /" [ 2pip)A{u) : A{u). 
JO Jn JO Jn 

The right-hand side members are estimated in the same way that for obtaining the estimate (j20p . We 
obtain 



Mp)|A(u)|2+ / / Vpip).V^ip) 



Jn 



<C(mo,po,S) + ^ r [ p|up + ^ r [ p\u + XJ\' + \A{n)\l f f pip). 

^ Jo Jn ^ JO Jn JO Jn 



(22) 



2.2.4 Gronwall argument 

Putting ([18]), ([201) and ([22]) together we get 



+ 2Q(p(r))) + r / M(rt|o(u)|" + 

/ JO Jn 



/i(p)|^(u)p 



+ 



T 

^n 



A(p)|div(u)|2+ 



C2M2 7o 



+ 



C0Q7 

2iV2 



/ 




Jn 





<C(mo,po,5) + C(|VS| 



Jn 



p|u|2+/ / p|u + U 



Jn 



Jn 



T 



+ ivui: 



Jn 



pip)- 



Since the gradient of both positive and negative powers of the density appear on the left hand side 



(recall that M 



< and iV 



7+n— 1 



> 0) and since the density is constant on d^l, we 
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can insure, thanks to Poincare, that J^ii{p) is controlled, via assumptions (jlip and (jl2p . by the 
pressure terms of the left hand side. Then, using a Gronwall argument, we deduce that all the left 
hand side terms of this last inequality are bounded and we can write the following estimates: 



|P lL°°(0,T;LT(f7)) 


< 


c(mo,po, 


u) 


L°°(0,T;L2(f7)) 


< 


c(mo,po, 


u) 


ll\/pVv9(p)||i«>(0,T;L2(n)) 


< 


c(mo,po, 


u) 


||iy^(p)Vu||i2((o,T)xQ) 


< 


c(mo,po, 


u) 


^ n — a — 1 , , 

2 ||L2(0,T;L2(n)) 


< 


c(mo,po, 


u) 


2 L2(o,T;L2(f7)) 


< 


c(mo,po, 


u) 



2.3 Approximate solutions and compactness 

The preceding a priori estimates are the key ingredient of the existence result. As soon as approximate 
solutions satisfy such estimates, compactness properties make it possible to extract a subsequence that 
converges to a weak solution of the initial model. The compactness arguments are exactly those given 
in the periodic case or the whole space in and more recently in the bounded case, see [8], that is 
why we will not detail it here. 

Let us just say some words about the sequences of suitably smooth approximate solutions to the 
compressible Navier-Stokes equations ([I|)-([2|) that preserve the estimates obtained in the previous 
section. 

The construction scheme of approximate solutions, using on additional regularizing effects such as 
capillarity, is provided in [6]. We thus introduce some modified Navier-Stokes equations for {pa,/3, Uq-,/?), 
always denoted (p, u) for sake of simplicity, depending on the regularizing parameters a and (3: 

dtp + div{pu) = 0, (23) 
dt{pu) +div{pu(g)u) - diva -rop\u\u- f3p\7{p'{p)A''p{p)) +aA^u = 0, (24) 

where the conditions ([7|)- (fT2]) are supposed to be satisfied. 

These regularizations allow to use some classical result in order to prove the existence of smooth 
solutions. The remaining work consists in showing that the additional terms depending on a and on (3 
lead to some weak solutions of our initial model ((l])-(l2]). 

Notice that we may not modify (|23p because the BD entropy is very closely related to the mass 
equation and some regularizing term in (123^ could cancel equation (llSh . For this model, since energy 
and BD identities are preserved, the stability arguments given in ^ and [7J lead to our existence result 
cited in the Theorem 11.21 

3 Sketch of proof of Theorem 11.41 (stationary case) 

The proof of Theorem 11.21 has been managed in the general case tq > 0. The only two points that 
have to be cleared in the stationary context concern the BD structure and the control of integral 
terms in the energy and BD formula. To control these terms we assume in this part, as announced 
in Theorem 11.41 that ro > 0. Moreover, in the three-dimensional case, we will assume the additive 
condition (fT3]) . 

These two additive conditions (the condition on tq and the condition (|13p ) will be used since in the 
stationary case we can not use Gronwall type arguments. 
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3.1 BD structure 



Bringing some modification to the mass equation could cancel the BD structure, that is why it is not 
clear that Subsection 12.11 can be directly adapted. For instance, it is dedicated to the failure for any 
semi-stationary model, whereas stationary conditions for both mass and momentum equation lead, 
following the same steps as for the equation (jlSp . to a similar equation: 

div(/)u (u + U)) + Vp{p) = div(2/i(/9)yl(u)) - rop|u|u, 
3.2 Control of integral terms 

In the stationary case, we will use the friction term to obtain a "good" estimate. More precisly in this 
case the energy estimate and the BD entropy estimate pT|) respectively write 



2^,{p)\D{vi)\' + / A(p)|div(u)|^+ro / p|u| 
n Jn Jn 



= - / (pu(g)u):Vu+ / 2fi{p)D{u) : D{u) +ro / /9|u|u • u, 
Jn Jn Jn 

[ 2Mp)|A(u)|2+ [ VP{p).V^{p)+ro [ p|u|3 
Jn Jn Jn 

= - /" (pu » (u + U)) : VS - [ 2ii{p)A(u) : A{u) 
Jn Jn 

+ ro / p|u|u-u — ro / /j|u|u-U. 

Jn Jn 



(25) 



(26) 



We estimate the terms of right-hand sides again (the nonhere detailed terms are exactly treated as 
in the nonstationary case). The constant C which appears does not depend on physical constants 
such ro, u... 



/ {pu u) : Vu 

Jn 



< ^ I p|u|3 + -|Vu|^ / p, 



4 Jn ro 



n 



• 


ro / p\u\u ■ u 






Jn 


- 4 Jn 



|3 
I oo 



n 



The only two terms which seem more difficult to control are the following 

Ti = / (pu0U):V5 and T2 = ro / p|u|u-U. 

Jn Jn 

• Using the definition of U and of 99, and using an integration by part, since div(u) = 0, we obtain 

Ti = / 2(u (g) V^(p)) : V5 = - / 2p{p) (Vu)"^ ■.Vu+ 2p{p) (u • V5) • n. 
Jn Jn Jan 

Since u • n = and u = u = on dO. we write 

(u • Vu) • n = Uj(9jUj)nj = Uj5j(ujnj) — Uj(5jnj)uj = — • Vn • u^ = I^u**), 
where II is the second fundamental form of dO.. By definition of ^(u) and -D(u) we have 

(Vu)^ : V5 = {D{u) - A{u)) : (Z)(5) + A(5)) = D{u) : D{u) - A{u) : A{u). 
Hence we obtain 

Ti< / p{p)\D{u)\' + \D{u)\l [ p{p)+ [ p{p)\A{u)f + \A{u)\l [ p{p) + 2p{pfU{u>^). 
Jn Jn Jn Jn 

10 



Finally, using the assumptions ([8]) and Q we can control fi with p as follows 

/ fi{p)= [ fl{p)lp<A+ [ p{p)lp>A<- [ p"lp<A + - / p"'lp>A<^^ + - / p"". 

Jn Jn Jn cq Jq c\ Jq cq ci Jq 

We obtain 

Ti< / p{p)\D{u)\'+ [ Mp)|^(u)P + 1^M1^ + ^ /■ ;,- + 2M/.)''II(u^). 
Jn Jn Co ci 

• For the term T2, since pU = 2\/p{p) and u • n = on 50, we obtain by integration by part 

72 = - / 2rop(p)(|u|div(u) + -^-(u-V)u) < 4ro / /i(p)|u||Vu|. 
V |u| / Jfl 

By the Young inequality, we obtain 

T2< f /i(p)|Vu|2 + 4r2 [ /i(p)|u|2. 
in Jn 

Using assumptions ([8]) and ([9]) , the fact that n > 2/3 and the Young inequality, we succesively deduce 
that 

4r2 [ /i(p)|u|2 < 4r2 [ p{p)\u\H + [ pip)\u\Hp^A 
Jn Jn Jn 



Co Jn Ci Jq 

4 ./o Cq 



Consequently we majore T2 as follows 



4 Jn eg cf Jn 



With the preceding estimates, the sum of the equalities ()25p and (j26p is written 
Mp)I^(u)P + ^ Mp)I^(u)|' + A(p)|div(u)|2 

With Cte = MlXl^ + illo^ + 2pipfU{n'). 

Co eg 

We conclude this section by showing that all the terms of right-hand side of the equation (|27p (except 
the constant Cte) can be controlled by the terms of the left-hand side. This result is due to the control 
of the density via the term p'^~^^'^\'V p\'^ ~ In\^ (P^) I where N = '^^^""^ . This term make it 
possible (using the Poincare inequality) to control p^ in H^{i}). ^From the Sobolev embeddings, we 
deduce a control of p in L'^^{i}) (for all q < +00 in the 2-dimensional case, and for all q < in the 
d-dimensional case, d> 2). 
If we assume that 

^m-2<qN (CI) 
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then all the integrals of the right-hand side of the equation (j27p are controlled (since m > 1, that is 
3m — 2 > m > 1). For instance, using the Young inequality, we write 

Ci Jn Jn 

where we can adapt the constant t^i such that the term 5i p'^^ is controlled by ^2N^ In 1^ (P^) f ■ 
3.3 Stability of weak solutions 

In the stationary case, the lack of estimates implies that the stability of weak solutions is conditionned 
by some specific profiles for viscosities and pressure. Some relations between the corresponding coef- 
ficients m, n, a and 7 may be considered. 

Let us consider a sequence of weak solutions Pk^u^ of the stationary equations ©-(HI). 
3.3.1 Estimates 

The preceding subsection leads to the following a priori estimates: 



llV/^(Pfe)Vufc 


L2(n) 


< 


c{Q, u). 


(28) 


||V'A(/Jfe)div(ufc) 


L2(n) 


< 


c(r2, li). 


(29) 


liv (e(p.)^O 


L2(Q) 


< 


c{Q, u). 


(30) 


l|V (pf ) 


LL2(n) 


< 


c{Q, u). 


(31) 


Wpk^l 




< 


c{Q, u), 


(32) 



where M = "i^-S^ < and iV = > 0. 

We are going to show that these estimates together with some compactness arguments lead to conclude 
that (yOfcjUfc) weakly converges to a solution (p, u) of the system (l3])-([l]). 

3.3.2 Compactnesses 

In order to cover the general case d G {2,3}, we will keep a coefHcient q such that i7^(r2) C L''(0) 
with continuous injection. In the d-dimensional case (with d > 2) we can choose any q such that 
q < 2d/{d — 2) whereas in the 2-dimensional case we can choose any q such that q < +00. In the 
sequel, we will denote by q such a real. 

• Compacity on the density - The estimate pip shows that the sequence p^ is bounded in H^{Q). 
Under the condition (|Cip and the fact that 3m — 2 > 1 for all m > 1, we have qN > 1. Consequently 
we obtain 

Pk^p in L^^'in). (33) 

In the same way, the estimate (j30p shows that the sequence p^^ is bounded in H^{Q,) (recall that by 
definition we have M < 0). We obtain 

— ^ - in L-'^^\Q). (34) 
Pk P 

We will note that —qM > 1 is satified in the 2-dimensional case taking q large enough and in the 
3-dimensional case taking q = 6 and using the assumptions given on page [3] for a and n. 
By the conditions ([8]) and ([9]), we obtain 

VpiPk) V pip) mL^(O), (35) 
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We will note that > 2 (using the condition (jCip and the fact that 3m — 2>m for all m > 1) and 
> 2 (since we previously prove that —qM > 1 and since n < 1). 



• Compacity on the velocity - On another hand, we know by (|28p that fi{pk)Vuk is bounded 
in L^(r2) and thus weakly converges in L'^{Q). /,From the identity 



Vufc = A//i(pfc)Vufc, 

we also conclude that Vu^ is bounded in U'{i}) with f = ^ — 2g]\7' note that r > 1 since we have 
previously proved that —qM > 1 > n. Moreover, since M < 0, we also have r < 2. 
Using Poincare inequality we obtain a bound for the sequence Ufc in W^'^'in). Thanks to the com- 
pactness W^'^'i^l) C L''{^}) for s < we obtain 

Ufc ^ u in L'{n), Vs < -^p^. (37) 

3.3.3 Limit 

In this subsection, we show that we can pass to the limit when k tend to +cxd for the nonlinear term 
in the equation ([3]) and (j4]). The "more nonlinear" terms in these equations are the following ones: 

Ti = P{pk), T2 = pk\uk\uk and Ta = div(^(pfc)Vufc). 

More precisely, the other nonlinear terms are div(/)fcUfc) and div(/>fcUfc(K)Ufc) which convergences (in the 
sense of distributions on $7) are consequences of the convergence of T2, and V(A(pfc)div(ufc)) which 
convergence is similar to the convergence of T3. 

• Convergence of the pressure term Ti - Recall (see assumption (|10p ) that the pressure is a sum 
of two pressures Ph+ Pc- 

■k Since Ph{pk) = o-pI convergence of Vph{pk) to Vph{p) in the sense of distributions on O 
comes from to convergence ([33]): 

qN 

pI^ p^ in L-{n). (38) 

We will note that ^ > 1. More precisely, this condition is satisfied in the 2-dimensional case 
(taking q large enough) and in the 3-dimensional case taking q = 6, and n and 7 satisfying the 
assumptions given pageO 

•k Then, we are interested in the convergence of the cold pressure term which writes as 

^Pc{pk)=Pk [Pk PcipkW Pk^{p^,<p*}j + Pk [Pk Pc{pk)ypk^{pk>p*}j- 

The only thing we have to obtain on the gradient of the cold pressure 'Vpdpk) is its boundedness 
in some L*(r2) space with t > 1. Recalling the assumptions (jlip and (I12p on the cold pressure, we 

know that /Ofc^"^"|Pc(PA:)Vpfc|l{pfc<p*} and /of "™''''^^''^Vc(Pfe)V/'fc|l{pfc>p*} are bounded in L'^{Q), 
respectively by ([30|) and ([3T]) . As a consequence, we can insure that Vp^Pk) is bounded in L*(J7) 
with t > 1 as soon as p^*^ ^ and p™^^'^^''^^ ^ are bounded in L'^{Q). 

Since — M — a < 0, using the convergence result we get the expected information "p^''^"" 
is bounded in L^(J7)" if we have 

2{a + M) < -qM. 
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We note that this condition is satisfied in the 2-dimensional case taking q large enough and in 
the 3-dimensional case taking q = 6. 

In the same way, since max{/3, 7} — > 0, using the convergence result ([33]) we get the expected 
information "p™'^^'^^''^^ ^ bounded in L^(r2)" if we have 

2(max{/3,7} - N) < qN. 

As previously, we note that the condition 2(7 — A^) < qN is satisfied as well in the 2-dimensional 
case as in the 3-dimensional case. Hence, we need the following condition 2(/3 — N) < qN which 
can be written 

4/3 < (g + 2)(n + 7- 1). (C2) 

2 1 

• Convergence of the friction term T2 - We write T2 = (/9fc|ufc|^)^ p^. Using the compacity (that 
is the strong convergence p3]) ) and the bound (see estimate ([32]) ) on /9fc|ufcp in L^{Q), we obtain 

2 3 

where / is the weak limit of (pfc|ufcp) ^ in L2 (J7). The last inclusion holds since qN > 1 (see condi- 
tion jnH)). 

To identify the limit /, we use the strong convergence for the density and the velocity: 
pI ^ inL— (17) and |ufcp ^ |up in L^(17), Vs < _ 

We deduce that / = /03|up if + < 1- We can show that this condition is satisfy in the 

two-dimensional case taking q large enough and in the 3-dimensional case taking g = 6 (and using the 
fact that iV > i, M < ^ and n < 1). 
Consequently, the friction term T2 satisfies 

21 2 1 
T2=Pfc|ufc|2= (pfe|ufc|3)^p3 ^p3|u|2p3 =p|u|2 inL^in). 

• Convergence of the viscous term T3 - Through ([28]) we obtain 



p{pk)Vuk g inL^(17). 



To identify the limit g, we use the strong convergence results das]) and ([37]). We get g = VMpJVu if 
we have the following condition 

m n , , 

This condition is satisfied in the 2-dimensional case taking q large enough. In the 3-dimensional case, 

^ and iV = ^ 

n 



taking g = 6 the condition (|39|) is written (recall that M = " ^ ^ and — n+j_\_- 



m < (7 + n — 1) ( 3 



1 + a — n 



Since n < 1 and a > 1, we have 3 — i^2-n ^ ^' deduce that the condition (j39p is contained in 
the condition (ICID in the 3-dimensional case. 



The viscous term /i(pfc)Vufc is written \/ fJ-{Pk){\/ p{Pk)^^k) which converges in L^{Q) if m < qN, 
condition which is a consequence of the condition (|Cip since 3m — 2 > m. We obtain 



div(/x(/9fc)Vuyfc) div(//(/9)Vu) in V'{Q). 
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3.4 Pressure and viscosity conditions 

Let's recapitulate all the conditions we need to get the integrabilities and compacities cited in the 
preceding subsections. Recall that (see assumption on page [3]) 

2 

7>1, a>l, m> I and - < n < 1. (40) 

The additional conditions (ICip and ()C2p write as 

3m-2<qN, 4/3 < (g + 2)(n + 7 - 1). 



In the two dimensional case, q can be chosen as large as we need, thus many inequalities are 
satisfied and Theorem 11.41 holds only with conditions (j40p . 

In the three dimensional case, q is any number smaller than 6. Thus, Theorem II .41 holds with the 
following additive conditions on the pressure and viscosities coefficients: 

1 

m<7 + n — -, (3 < 2(7 + n — 1). 
These conditions exactly correspond to the condition (jl3p . 



4 Behaviour in thin domains 

In this part, we derive the compressible Reynolds equation. Formally, this equation comes from the 
compressible Navier-Stokes equation in a thin domain, that is when one of the length is assumed to 
be smaller than the other directions. The main applications of this kind of behavior relate to the field 
of lubrication (see the Introduction). Within such a framework, the thin domain is of the following 
form 

= {(x, z) G M.'^-^ X M ; x G O C M"'"^ and 0< z < eh{x)}, 

where O is a bounded domain in M'^~^ and the height : O ^ M is a regular and periodic function. 
Note that to be able to define a periodical function, the domain O must be rectangular. In the case 
of the dimension d = 2 this is not a resctriction. In the case of the upper dimension, this situation 
corresponds to realistic physical situations. Moreover, it is possible to consider other conditions on the 
lateral boundaries. For all these aspects, consult thesis of S. Martin [T7], as well as named references. 
We assume that h > /imin > and up to a normalization, we can assume that /imin = 1- The size of the 
bounded domain O C W^^^ is assumed to be of order 1. The non-dimensional number e corresponds 
to the characteristic ratio between the characteristic lenghts of O and the characteristic height eh. 




u = (y, 0) p = pb X G R'^- 



Boundary conditions on 0,^: According to the results of the preceding parts (and according to the 
periodic results, see for instance ^), the boundary conditions which we impose are the following 
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(i) - Periodic conditions for the velocity and the density on the lateral boundaries (i.e. for x G dO). 

(ii) - Dirichlet conditions for the velocity on the top and bottom surface 

u = {V, 0) G M'^-^ X M for z = 0, u = (0, 0) for z = h{x). 

(iii) - Constant density on each connex component 

p = Pb for z = 0, p = pt for z = /i(x). 

The goal of this part is thus to justify in a rigorous way the compressible Reynolds equations, i.e. to 
determine the limit when e tends to of the stationary compressible Navier-Stokes equations ©-([I]). 

Remark. The result which is shown here concerns the justification of the Reynolds equation from the 
stationary Navier-Stokes compressible equations. Of course, the same method would allow to give a 
rigorous justification of the non-stationary Reynolds equation from the non-stationary Navier-Stokes 
equations. 

4.1 Rescaled equations 

In such a domain the unknowns of the equations ©-(HI), i.e. velocity and density, depending on e 
are denoted with a subscript and ps- The first stage consists in rewriting these equations ([3])-(|4]) 
in a domain independent of e. For that, we introduce the change of variable Z = z/e. We define the 
rescaled domain 

n = {(x, Z) £ M'^-^ X M ; X € O C M'^"^ and < Z < /i(x)}. 

In worries of simplifications, computations and notations used later will be made in dimension d = 2. 
In this case, the impose velocity ^ is a real number, which will be assumed to be positive: ^ > 0. If 
the three dimensional case {d = 3) is really different (for instance when we use the classical Sobolev 
injections) we shall apparently refer to it. 

In the studied context (for example that of lubrication), we know that the pressure depends on the 
thickness e of the domain as 1 /e^ (see |17) and the cited references) . We define a normalized pressure 
by Pe = e^Pe- In the same way, if horizontal velocity is of order 1 (this order of magnitude depends 
in fact on the size of the velocity imposed on the boundaries of the domain, here we suppose that V 
is of order 1) then vertical velocity will be of order e: = {ve^ewe). We can obtain the following 
equations in Vt: 

d,{peVe) + dz{peWe) = 0, (41) 
PsVed:r.Ve + PeWedzVe = 2dx{p{pe)dxVe) + -^dz{p{Pe)dzVe) + dz{p{Pe)dxWe) (42) 

+dx{KPe){dxVe + dzWe)) - ^dxP{pe) - r^Peivl + e^w'^^)^Ve, 

1 2 

PeVedx{eWe) + PeWedz{eWe) = -dx{p{p£)dzVe) + edx{p{pe)dxWe) + -dz{p{Pe)dzWe) (43) 

11 1 

+ -dz{\{pe){dxVe + dzWe)) - —dzP{pe) - £roPe{vl + E^wlf^We- 
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As we have seen it in the Subsection 12.11 we are also interested in other particular forms of these 
equations. Refering to (jl5p and noting Ue = (V^jVF^) = (^I^^^iEil^ ^^^^^^^-^ ), we can write 

PeVsdx{Vs + Ve) + PsWedz{Ve + = dz{p{Pe)dxWe) " \dz{p{Pe)dzVe) 

-^dxP{p,) - roPe{vl + e'wD'^Ve, (44) 

PsVsdx{eWe + We) + PeWedzieWe + Ws) = edx{p{p£)dxWe) - ^dx{p{pe)dzVs) 

-^dzP{pe) - erapeivl + e^wl)^We. (45) 

4.2 Convergence of the compressible Navier-Stokes equations to the compressible 
Reynolds equations 

Theorem 4.1 A solution {pe,Vs,£W£) of system li41\ )- (4^ with conditions ([7D- (fi^[ J satisfying the 
preceding boundary conditions (i), (ii) and (Hi) converges to {p,v,w) inU^{Q) X {U^{Q)) ,for some 
ri,r2 > 1, when £ tends to 0. 

At the limit, the following system holds in D'(n); 

nh 

dx{ / pvdZ) = 0, (46) 
Jo 

-dz{p{p)dzv) + dxP{p) = 0, (47) 
dzP{p) = 0. (48) 

Moreover, the horizontal velocity v satisfies the boundary conditions v\z=o = V, v\z=h(x) — '^'^^ ^ 
and p are periodic with respect to the x variable. 

Remark. 

1. The boundary conditions on the density are not conserved through the limit e — > and this is 
essential to get a non constant density p in the limit system. 

2. The limit model is not coupled with the vertical velocity any more. However, we have a limit 
equality for the weak limit w of ewe'- dz{pw) = 0. 

3. Notice that if P' is not zero almost everywhere then dzP{p) = <^=^ dzp = 0. As in the 
incompressible case, we can integrate twice the first equation ()47|) with respect to variable Z 
and use the condition (jl6]) in order to obtain: 

4. It is important to notice that for such a Reynolds equation, maximum principle is proved (see 
for instance pS]). Consequently, there exists a constant pram > such that the solution p of 
the Reynolds equation ([Mj) satisfies p > Pmin- We deduce that for e small enough, we have 
Pe ^ > 0) that imply that assumptions ^ and ([TT]) are useless. 

5. The proof presented here can easily be extended to the nonstationary case (by using the result 
of the Theorem ll.2p . We thus justify the nonstationary compressible Reynolds equation 

as the limit of the compressible Navier-Stokes equations in thin domain. 
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5 Proof of Theorem 14.11 



For all £ > 0, the existence of a suitable solution of (j4ip - (j43p is given by Theorem 11.41 coupled with 
classical results in periodical cases (see [7]). So, let us consider (^gjUg) such a solution. 
The aim is to obtain estimates of {pe, u^) which are on the one hand non-depending on e, on the other 
hand sufficient to pass to the limit e — > in the equations (|4ip - ()43p . The main difficulty comes from 
the non-linearities which require strong convergences of some terms. 



5.1 A priori estimates 

To write energy estimates, we take again the method of previous sections. Recall that this method 
requires the introduction of a velocity lift u^. Within the framework which interests us here (see 
Figure page [T5|) . the velocity lift that we will use is: = ('Ue,!!;^) with 



V{1- Z) if < Z < 1 
if Z > 1 



and Wi; = 0. 



Note that it is very easy to regularize this velocity field, keeping the form = (vs{Z),0). 
Lemma 5.1 The energy and BD formula write as follows 



piPe)\d^Ve\'^ + 2 / p{pe)\dzW,\'^ + / pipe) 



edxWe + -dzVe 
£ 



+ / \{pe)\d^Ve + dzWe\''+ro [ Pe{v^, + £W,f/^ = Sf 

Jn Jn 



I 

Jn 



edxWe 



-dzVe 



+ 



1 Co 

C2M2 
1 ^007 



M\2 



+ 



1 Co 



£4 C2M2 



\dz{aps)r\' 



where 



Sf = Ti + - [ p{pe) i- dzVe + £ d^We) d,Ve + ro [ Pe Ve Ve {vj + £^wlfl^, 

£ Jn ' Jn 



with Ti= PeVedxVeVe+ / PeWedzVeVe, 



S^ = Sl + 2 / d^p{pe) We dzVe + 2 Tq / [ve S^nipe) + We dzfJ-iPe)) {vl + £ 

Jn Jn 



'wlYi\ 



(50) 



(51) 



Proof of the energy estimate (I50|) 

We first note that the function permits to get homogeneous Dirichlet boundary conditions on the 
corrected velocity — Ug . To get a first energy identity, we sum equation (j42p multiplied by the new 
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(52) 



horizontal velocity Ve — and equation ([33|1 multiplied by e{w£ — We) = eWe- We obtain 

/ PeVedxVeiVe - Vs) + / PeWsdzVe{Ve - + e'^ / PeVeWed^We + / Pew'^dzWe 

Jn Jn Jn Jn 

= 2 / d:^{p{pi;)d:,;Ve){Ve -Vi;) + ^ [ dz{p{Pe)dzVe) {Ve - Ve) + [ dz{p{Pe)dxWe) {Ve - V^) 

Jn £ Jn Jn 

+ / d^{p{pe)dzVe)we+ / 9^ (^(/Je)^^;^^) + 2 / dz{p{Pe)dzWe) 

Jn Jn Jn 

+ / d^[\{pe)(d^Ve + dzWe)){Ve-Ve)+ / 5^ (A(/9£) ((9^.U£ + ^zWe)) 

Jn 

- \ I d^P{pe){Vs - Ve) - \ I dzP{Pe)We 

Jn Jn 

-ro Pe\v'^ + e'^w'^l^VeiVe - Ve) - Tqe'^ / Pe\v'^ + e'^w'^\^ w'^ . 

Jn Jn 

By the complexity of this equation, we will deal with these terms by group. The terms of the first 
line (left hand side of the equation) are known as convection terms. Those of the three following lines 
will be called viscous terms. The fifth line corresponds to the terms of pressure whereas the last line 
contains the friction ones. 

Convection terms - We first use integrations by parts (without boundary terms thanks to the 
x-periodicity and the boundary conditions for We) and the divergence free conditions: 

/ PeVedxVeVe+ / PeWedzVe Ve = - / PeVe Ved^Ve - / PeWeVedzVe- 

Jn Jn Jn Jn 

This contribution is zero (since it equals to its opposite). In the same way we have 



e / PeVeWedxVUe + e / PeWedzWe = 0. 

Jn Jn 

Finally, all the terms of the left hand side of (I52|) disappear except those containing Ve, denoted Ti: 



Ti = / PeVedxVeVe+ I PeWedzVeVe- 

Jn Jn 

Viscous terms - These terms are easily computed using integrations by parts. In any integration 
by parts, no boundary integral term appear thanks to the vertical correction induced by and the 
periodicity in x. 

Pressure terms - We also remark that, using div{peUe) = dx{peVe) + dz{peWe) = and dxVe = 0, 
the pressure contributions vanish. In fact, noting li-^Pe) = ^—y^-, we have 

f dxP{pe){Ve -Ve) + [ dzP{Pe)We = - I Ii{pe){dx{peVe) + dz{PeWe)) - [ P{pe)dxVe = 0. 

Jn Jn Jn Jn 

Friction terms - Clearly , the friction terms (that is the terms containing the friction coefficient tq) 
appearing in estimate (j52|) write 



-ro / Pe\Ve + e'^w'^l^ +ro / Pe\v'^ + e'^w'^l^^VeVe- 
Jn Jn 



All these calculations give the first identity of Lemma 15.1 
Proof of the energy estimate (15 ip 
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(53) 



Refering to Section [21 we can also write another equality related to the BD entropy. To get this BD 
formula, we sum equation (HH) multiplied hy Vgr — Vir + Ve and equation (|15]) multiplied by We — 'W£ + We 
(for information, recall that Ve = and We = ^^^Mpd). We obtain 

/ PeVedxiVe + Ve){Ve -Ve + Ve) + PeWedz{Ve + Ve){Ve - Ve + Ve) 

Jn Jn 

+ / PeVedx{eWe + We){eWe + We) + / PeWedz{eWe + We){eWe + We) 

Jn Jn 

= / dz{p{Pe)dxWe){Ve -Ve + Ve) 3 / dz{p{p£)dzVe){Ve -Ve + Ve) 

Jn £ Jn 

-e I dx{p{pe)dxWe){eWe +We) + - [ dx{p{pe)dzVe){eWe + We) 

Jn £ Jn 

- \ [ dxP{pe){Ve -Ve + Ve)- \ [ dzP{Pe){eWe + We) 

Jn e"* Jn 

-ro pe\v^ + e'^W^\^Ve{Ve - Ve + Ve) - ro / epe\v^ + e'^UJ^\^ UJe(eWe + We) . 

Jn Jn 

In addition to those which are common with the energy, we deal with every terms: 

Convection terms - Let us look at the first line of the equation ()53p . The terms which are not 

already treated to obtain identity (jSOp are the following ones: 

/ PeVedxVeVe + / Pe Ve dxVe {Ve - Ve) + / PeVedxVeVe 

Jn Jn Jn 

+ / PeWedzVeVe + / Pe We dzVe {Ve - Ve) + / PeWedzVeVe- 

Jn Jn Jn 

The divergence free condition dx{peVe) + dziPeWe) = is also strongly used here. With the boundary 
conditions (periodicity with respect to x and boundary conditions on tt; for Z G {0, 1}) we deduce 

that all terms are equal to zero by integration by parts, except this one: / peWeVe dzVe- 

Jn 

Since peVe = 2dxp{pe), we write this additive term as 2 / dxp{pe) We dzVe- 

Jn 

Viscous terms - The only additive terms compared the first energy identity (see equation ()50p ) are 
the following 

/ dz{p-{Pe)dxWe)Ve - \ [ dz{p{pe)dzVe)Ve - E [ dx{p{pe)dxWe)We + - [ dxipipe)dzVe)We 

Jn ^ Jn Jn £ Jn 

= / p{pe)dxWe{edxWe - dzVe) + \ fl{pe)dzVe{dzVe - edxWe) = 0. 

Jn ' ^ ' Jn ' ' 

=0 =0 

The equality edxWe — dzVe = comes from to the following computation 

edxWe-dzVe = 2dxi^-^^)-2dz(^^^)=2dzp{pe)dx(-)-2dxp{pe)dz(- 
^ Pe ^ ^ Pe ^ ^Pe^ ^ Pe 



-) - dxPedzi- 

Pe' ^Pe 



2p'{pe) dzPe dx[ — ]- dxPe dz[ — ) = 2 ^^P^ ^''P^ ~ ^^P<' ^^P^ 

Pe 



2m'(p.) 



0. 



Pressure terms - For the pressure terms, we can rewrite what we wrote for the energy, say, the tests 
against u^ — u^ are equal to zero thanks to (I4ip and dxVe = 0. The remaining part of the pressure 
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contributions have been discussed and expressed in (|T8j) . let's recall it here: 

\ [ d,P{pe)Ve + 4 / dzP{Ps)We 

Jn Jn 

s'' Jn ^ ' Pe £ Jn ^ ' Pe 



1 Coa7 [ ^ 2 , 1 Coa-f f 



Friction terms - Let's now deal with the friction terms, in addition to those which are common with 
the energy, we have to say some words about the ones which are specific for the BD formula. Using 
the definition of and W^, these additive terms are written 

ro / Pe {V^e + W^f^^ {Ve Vs+SWe W^) = 2 Tq / {v^, + wl^'^ [v^ d^Pe) + We dzp{pe)) ■ 

Jn Jn 
Putting all these inequalities together, we obtain ([5T|) . □ 



The left hand sides of the estimates ()50p and (j5ip make appear the L^(J7)-norm of the following terms 

PeldxVel"^, Pe\dzWe\'^, e"^ PeldxWel"^ , ^\dzVe\^ , \e\dxVe\^ , \e\dzWe\^ , rQPe\Ve\^ , e^rQPe\Wef., 



e'pe\dxwe\\ ^\dzve\\ ^,\dx{p^)? and 

where we have noted pe = P'iPe)-: = ^{Pe) and = ^{pe)- We now prove that the right hand sides 
of the estimates (j50|) and ([5T]) . that is the quantities Sf and S*!, can be controlled by such terms. 

Control of 5^^ - /,From the definition of Sf, we express Sf as follows: Sf = Ti +T2 + T3. 
The contribution Ti is writterU 

Ti= Pe Ve dxVe Ve + Pe We OzVe Ve 



{pl/'ve) {pT'^dxVe) {pf-'-^"ve) + / {pl'^SWe) (^^^) (p(4-3™)/6,~ ) 



n 



+ '-^1 Pewl + ^J p7\dzVe\' + 4r, 

4 Jn e Jn 



\"e\oo I Pe ' 



In this last inequality, 5 can be choosen as small as possible (due to the Young inequality). Moreover 
the constant C does not depend on the physical constants e, tq, ci... nor 5. The term 5 J^p'^\dxVe\^ 
is bounded using the assumptions ([9]) on /x as follows: 

5 / pT\dxVe?=6 I p^\dxVe\^\p-^A + 6 [ pT\dxVeflp<A 

Jn Jn Jn 

< —J^P{pe)\dxVe\'^tp:^A + ^^;^ J^p';^\dxVe\'^lp<A 
<— [ pipe) IdxVel"^ lp>A + [ pipe) IdxVel"^ tp<A 

ci Jn Co A Jq 



^Note that this term is not treated as the corresponding one in the proof of Theorem ll.4l In the proof of Theorem II. 4 1 
we write Ti — — dzVe- This is correct but inappropriate in our situation. In fact, we want to obtain estimates 

with respect to the parameter e and we have no control on We, but only on eWe and on \dzVe. 
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Taking 6 = minj^, ^o^™ " | obtain 



6 I pT\d.Ve\'<l [ KPe)\d.V,\\ 

n Jn 



In the same way, we estime p^ldzVsl"^ and we obtain (note also that |fe|oo is bounded by 1) 

C 



.2p LPs 



4-3m 



Since \dzVe\ < V , the contribution T2 is controhed as fohows 

1 /■ 1 _ 

T2 = - p{pe) (- + £ d^We) d^Ve 

^ Jn £ 
F /■ 1 

< — / ;U(/9e)|-(9^t;£: +e5a;U;£| 

< ^ / p{Pe)\^ dzVs + £ d^We\'^ + ^ j p{pe)- 

Using assumptions ([8]) and dH), we deduce that 

If ,1 ,2 y2 /■ 

7^2<:t / p{pe)\-dzVe + ed,W,\^ + —^ + —^ / Z^^- 

2 Jf^ ' e 2co 2ci 

The term is treated using the Young inegahty (in the fohowing form: AB^ < 6AB^ + ^ys^A, for all 
6 > 0). We obtain 

T3 =ro / p,^ye^^e(t''+'^'^«?)'^' 

Jn 



<ro I Pe\Ve\^+ro I Pe\Ve\\£We\ 



<^ Pe \Vs\ + — I ps \£We 



n 

2 



<? /' Pe\Vef+"-^ [ Ps\w,f + Cro [ p 

Jn Jn 



n 



t3 

Here, the constant C does not depend on e nor on tq. 

Control of 51 - ^From the definition of Sf, we write Sf as follows: S*! = S'f + T4 + T5. 
Using an integration by part and {dzVel < V, the contribution T4 is controlled by 

-2 



T4 < 2V j^p{pe)\d^We\ <^J Pipe) + y J^PiPe)\d.W. 
Using assumptions ([8]) and we deduce that 

T4 < - / p{peMWe\^ + 2 + 2 / 

2 CO ci Jn 



2 



Finally, the term T5 is written 

n = 2ro [ {Ve d^p{pe) + We dzp{pe)) {vl + £^ wlf'^. 

Jn 

For sake of simplicity, we only treat one example of this contribution (more precisely, the term 
2ro f^v'^ dxp{pe)), the other terms are similar. We have 



2ro / vl dxp{pe) = -^r^ / p{Pe)VedxVe / p{pe)\dxVe\^ +^rl / p{pe)\v.,^'^ 



n 

< \ / P{pe)\dxVe? +^rl I p{Pe)\Ve?'^p<A + ^rl / p{p,)\Ve\' 1 p>A- 
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Using succesively the assumption ([8]), the fact that n> 2/'i and the Young inequahty we obtain 



Jn Co Jq Co Jq 4 Jq 



In the same way, we use the assumption ([9]) and the Young inequality to obtain 

r'o [ PiPs) \Vs\' 1,>A < ^ / \V,\' <"i I Pe\Ve? + ^ f pT~' ■ 
Jn Cl Jq 4 Jq cf Jq 

Here again, the constant C does not depend on e, ro, co or ci. 
Finally, the term is controlled by terms like 

T,<\l p{Pe)\d.V,\^+"i I Pe\Ve? + ^ [ pg^^ + ^ ^0 f . 

Jn ^ Jn Cl Jn Cq 

Estimates - The result of Lemma 15.11 and the control of the terms Sf and allow to obtain the 
following estimates (note that in this estimate and in the following ones, the only constant that we 
will reveal will be e, the others will be taken equal to 1 to simplify calculations). 

/ p{ps) \dxVe\'^ + I p{pe) \dzWe\^ + p{pe) + 4 / p{Ps) \dzVe\^ + I Pe \Ve? 

Jn Jn Jn ^ Jn Jn 

+ e' I PeW? + y, [ |a.(ef)P + ^ / |5z(Cf)l^ + ^ / mP^)n^, [ \dziP^)f (54) 
Jn ^ Jn e Jn e Jn e Jn 

< I pt''"' + \ + ^ [ pT+ I Pe+ I pf'-' + l. 

Jn ^ ^ Jn Jn Jn 



The right hand side of the estimate (j54p make appear different powers of p^- We must show that all 
these terms can be absorbed by some left hand side terms, taking e small enough. 
Since m > 1, we have 3m — 2 > m > 1 and 3m — 2 > 4 — 3m. So, we must control the term p\^~'^ 
and the term when 4 — 3m < 0. 

• With the Poincare inequality, a control on |V(/9^)|^2(q) implies a control on Due to the 
Sobolev embeddings, this allows to control \p^\Li{n) for all q > -j^. 

Consequently, if 

3m - 2 < giV (Cl) 

then the term Pe™'"^ can be absorbed by the terms ^ l^zip^)]"^ and ^ l'9x(p^)P as soon as e 
is small enouglu. 

• Using the same arguments, from the definition of the function ^, the term \ dz{C^)\'^ allows to 
control JqPs'^, and recalling that M < we deduce a control of l/p^ in L~'^^{Q). Consequently, 
when A — 3m < we can absorbed the term Pe~^^ by the term \ l^^(Ce^)P as soon as e is small 
enough, under the condition 

qM <4- 3m. (C3) 
Finally, under the assumptions 3m — 2 < qN and qM < 4 — 3m we obtain for e small enough: 

/ P{Pe) lOxVel"^ + I p{Pe) \dzWe\^ + ^ / ^^iPe) \dxWe\^ + ^ / ^^iPe) \dzVe\^ + [ Pe \Ve\^ 

Jn Jn Jn ^ Jn Jn 

+ e' [ Pe\Wsf + -, [ \dx{aPerr + -, [ \dz{aPer)\' + -, [ Idxip^^ + \ [ |5z(pf)P<^. 

Jn £ Jn ^ Jn ^ Jn ^ Jn ^ 

''if we assume that 3m — 2 < qN, that is the condition HC1|I . then it is possible to control p^"""^ without taking e 
small but just using a Young inequality, see for instance the proof of Theorem [L4] where such a method is used. 
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We deduce the following bounds: 



\\\/ Ll(pi;)dxVe 


Il2(S1) 


< 


1/e 


(55) 


II V Ll{pp)dzVp 


T 2 ('O'l 


< 


1 


(56) 


\\-\/ IJL{Pe)dxWe 


L2(r2) 


< 


l/e' 


(57) 


\\\/ ^J'{pe)^zWe 


L2(f7) 


< 


l/e 


(58) 






< 


1 


(59) 






< 


e 


(60) 




LL2(n) 


< 


1 


(61) 






< 


£ 


(62) 


3 


L2(f7) 


< 


l/e 


(63) 


3 


L2(r2) 


< 


, 5 
1/62 


(64) 



where M 



n— a— 1 



< and N 



5.2 Compactness on the density 

Let us recall that, for simplicity of the notations, the preceding calculations were carried out on a 
two dimensional domain, i.e. on J7 C M"^ with d = 2. Of course, the latter remains valid in higher 
dimension, in particular in dimension d = 3. In this paragraph, we will strongly use Sobolev injections 
which are depending on the dimension. In order to cover the general case d S {2,3}, we will note q 
any real such that H^{i^) C L'^{^) with continuous injection. 

The previous estimates show that pf^ and p^ are bounded in H^{Q,). Therefore, we can write that i^ 
qN >1 and -qM > 1 then 

Pe^p inL«^(0), (65) 

p-i^p-^ inL-«^^(17). (66) 

Refering to the conditions ([8]) and ([9]), we also conclude that, for all q also satisfying 2qN > m and 
-2qM > n 

(67) 



VKPs) Vpip) in L - (O), 



VPiPe) 

5.3 Compactness on the velocity 



-2qM 

in L " (17). 



We know by ([56|) that yp(j)e]dzVe is bounded in L?'{^1) and thus weakly converge in L^(J7) to some /. 
^From the identity 

dzVe = ^—— y^p{pe)dzVe, 
VPiPe) 

we also conclude that dzVe is bounded in L''"(J7) with ^ = ^ — 2p7- 



^Under the condition (|C1|) and the fact that 3m — 2 > 1 for all m > 1, in fact we always have qN > 1. In the same 
way we will note that —qAI > 1 is satified in the 2-dimensional case taking q large enough and in the 3-dimensional case 
taking q = 6 and using the assumptions given on page |3] for a and n. 



As previously, these two conditions are satisfied, using the condition (|C1|) . the fact that 3m — 2 > m for all m > 1 
and "^"^"""^ > 2 (since we previously prove that —qM > 1 and since n < 1). 



24 



We note that r > 1 since we have previously proved that —qM > 1 > n. We have 

dzVe dzv in V^{Vt). 
Remark that we necessarily have r < 2. 

As for the derivatives with respect to Z, we know, using the bound (j55p and those which come from 
the convergence (I68p . that edxVs is bounded in U'[^) and thus weakly converges. Thus, using the 
Poincare inequality, we also get the bound of in U'{0,) and then 

vs^vin (70) 

ev,^OmWl''{Q). (71) 

Taking the same way, we also have 

eWe -^win L;(0), (72) 

edzWe dzw in (73) 

e^We ^0 inW^'^'in). (74) 
Moreover, thanks to the compactness W^''^{Q) C L'^{Q) for all s < r' where = f ~ 3 we can write 

£7;^ ^ in L"(J7), Vs < (75) 

a — r 

s^We^OmL'{n),ys<^P^. (76) 
5.4 Limit in the momentum equation 

Let us rewrite the two components ()4ip and (j43p of the momentum equation: 



e PeVed^Ve + £ PeWedzVe = 2e dx{fl{pe)d:^V^) + dz{fJ'{Pe)dzVe) + E dz{li{Pe)dxWe) 

+e'^dx{X{pe){dxVe + dzwe)) -d^{P{pe)) +rQe'^ Pe{vl + e^wl)^ve, 
e^PeVedxWe + e^p^WedzW^ = dx{f^{ps)dzVe) + e^dx{p{pe)dxWe) + 2e^dz{f^{p£)dzWi;) 

+e^dz{Xipe){dxVs + dzWe)) -dz{P{Pe)) +rQe^ Pe{vl + e^wD'^We. 

We initially will show that the bold terms admit limits when e tends to 0, then that all the other 
terms tend to zero. The method is exactly the same one as that developed in part 13. 3. 31 We just will 
specify the dependences in the parameter e. 

• For instance, putting together assumptions on the pressure (fTT|) - (fT2]) . the estimates (j59P - (j62]) and 
the strong convergences of the density (f65]l - (f66l) . one obtains the convergence of dx{P{pe)) to dx{P{p)) 
and the convergence of dz{P{pe)) to dz{P{p)) in the sense of distributions on $7, since we introduce the 
same hypotheses on the coefficients, see for instance the convergence of the term Ti in subsection 13.3.31 
and the assumption (|C2p . 



As for the term T3 in the subsection 13.3.31 under the condition (139]) we have 



^/p{Pe)dzVe ^ ^/p{p)dzv m L^Q) . (77) 

Since m < qN, we end to the convergence: 

p{pe)dzve ^ p{p)dzv in P'(17). (78) 

We are now going to show that the other terms tend to in the sense of distributions when e goes 
to 0. 
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• As for (j77j) and ([78|) . we also get, for the x-derivatives, 

e^^L{pe)dxVe ^ inL2(17), (79) 

e^(p,)a^t;, ^ inP'(17). (80) 

Therefore, we obtain the convergence of 2e^ dx{P'{pe)dxVe) ^ in P'(r2) and also dx{X{pe)dxVe) —>■ 
and e'^dz{)^iPe)dxVe) since A(s) = 2{sp'{s) — p{s)) satisfies the same integrabilities as p. Taking 
the same strategy as for the convergences already obtained through (f77|) . (fTH]) . (f79]) and (fHOl) . we finally 
can give the ones of the vertical velocity : 



e^J p{pe)dzwe p{p)dzw in L^(ri), (81) 

ep{pe)dzwe^ p{p)dzw inP'(O), (82) 
e^^yp{pe)dxWs ^ inL2(17), (83) 

eV(/o.)5xW^. ^ inV'in), (84) 

which answer the questions of convergence for the viscous terms containing Ws- 

• For the friction terms, we will just do it for the term tqe^p^vI (the other terms may be treated using 
exactly the same way, moreover the convective terms are also treated in the same way). Using the 
same method as in Subsection 13.3.31 in particular when we have treated the term T2, we have 



2 2 -/2 Sn-i 

roe PeV^ = roe 3 (e peV^)3 pf 



1 



^Prom the strong convergence of p^ and the weak convergence of e^peV^ we deduce that roe^p^vl tends 
to 0. Thus, we can insure that (j4ip -(|42 p converge to ()46p ~ (j47p in the sense of distributions. 

5.5 Limit in the mass equation 

To pass to the limit e — > in the mass equation dx{peVe) + dziPeWg) = 0, the main difficulty comes 
from the fact that the vertical velocity Ws does not have a limit. We thus will use the following 
equivalent form 

dx( PsVsdz) = 0. 



^ Jo ' 

It is an equivalent form in the following sens: if a velocity such that dxi^ J PeV^ dZ^ = exists 
then we can build a vertical velocity Ws such that dx{peVe) + dz[peWe) = and ^£1^=0 = We\z=h = 0. 



1 /"^ 

This is enough to define u)e = — / dxipeVe 

P h 

To prove that pi^v,, tends to pv we write 



1 ^ 

PeVe = {pevlY pi ■ 



Using the strong convergence of pe and the bound on Pkv\ in L^(0), we obtain (as soon as qN > 1, 
which is implied by the condition (jCip ) 

PeVe^gp^ inL^(ri), 



1 



where g is the weak limit of (/OeIue^) ^ in L'^{^1). To identify the limit g, we use the strong convergence 
for the density and the weak convergence for the velocity: 

pl^pi inL3«^(J^) and ^ v in L'^(O). 

We deduce that g = psv \i ^ + 3^ < 1, condition which is implied by the condition ([39]) . 
We deduce the following convergence 



dx ( ^ PeVe dZ^ ^ ( ^ dZ^ in ^'(^) • 
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5.6 Pressure and viscosity conditions 

In addition to the conditions given in Subsection 13.41 we have to suppose that (see assumption ()C3p i 

qM < 4 - 3m. 

In the two dimensional case, q can be chosen as large as we need, thus many inequahties are 
satisfied and Theorem 14.11 holds only with conditions (j40p . 

In the three dimensional case, q is any number smaller than 6. The additive conditions on the 
pressure and viscosities coefficients can be summarized by: 

7 

m < a — n -\ — . 

3 

This condition exactly corresponds to the condition (jl4p . 
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